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Abstract 

In this note, we give a closed formula for the partition function of the dimer model living on 
a 2 x n strip of squares or hexagons on the torus for arbitrary even n. The result is derived in 
two ways, by using a Potts model like description for the dimers, and via a recursion relation 
that was obtained from a map to a ID monomer-dimer system. 

The problem of finding the number of perfect matchings can also be translated to the 
problem of finding a minmal feedback arc set on the dual graph. 



1 Introduction 



In this note, we give a closed formula for the partition function of the dimer model living on 
a 2 x n strip of squares or hexagons on the torus for arbitrary even n. 

The dimer model is concerned with the statistical mechanics of close packed dimer ar- 
rangements on a bipartite graph. The real-world representation of the dimer model is the 
adsorption of diatomic molecules on a crystal surface. 

In the 1960s, the question of how many perfect matchings exist on a plane graph was 
solved independently by Kasteleyn [1, 2], and Temperley and Fisher (3HH: the total number 
is given by the Pfaffian of a signed, weighted adjacency matrix of the graph (the Kasteleyn 
matrix). Much of the original interest in the dimer model arose because it provides a simple 
and elegant solution for the 2-dimensional Ising model EQ. 

The problem of enumerating perfect matchings is of course a classical problem in graph 
theory and combinatorics (see e.g. [6 |), and can also be phrased in terms of domino tilings 0. 
During the last years, the interest in the dimer model was revived thanks to its manifold 
connections to other branches of mathematics and physics, such as the topological string 
A-model 13191/ real algebraic geometry |H0lHT] , BPS black holes from D-branes wrapping 
collapsed cycles [12] and supersymmetric quantum mechanics and categorification tech- 
niques Ifl3| . Furthermore, a correspondence between the dimer model and quiver gauge 
theories arising from D3-branes probing a singular toric surface was discovered and worked 
out in great detail lT4l[T5l[T6l[T7l[T8| . An explanation of this correspondence via mirror sym- 
metry was given in [ 19 1. 

The plan of this note is as follows. We briefly introduce the dimer model and give a 
Potts-like description of the dimer model living on a 2 x n strip on the torus. Using this 
description, we derive a closed formula for the Newton polynomial for any value of n. The 
same result can also be derived with a recursion relation obtained by mapping the problem 
to a one-dimensional monomer-dimer system, and is given both for a strip of squares and a 
strip of hexagons. 

Furthermore, the question is translated to the problem of finding a minimal feedback arc 
set on the dual graph. 

A bipartite graph ^ is a graph in which all vertices can be colored black or white, such 
that each black vertex is only connected by links to white vertices and vice versa. Let M be 
a subset of the set E of edges of . M is called a matching, if its elements are links and no 
two of them are adjacent. If every vertex of & is saturated under M, the matching is called 
perfect. Such a link that joins a black and a white vertex is called a dimer. The dimer model 
describes the statistical mechanics of a system of random perfect matchings. In the simplest 
case, we ask for the number of close packed dimer configurations, i.e. the number of perfect 
matchings. 
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Kasteleyn HJI3 introduced an orientation on , which leads to a signed adjacency matrix 
K, now called the Kasteleyn matrix. The Pfaffian of K gives the number of perfect matchings. 
A Kasteleyn orientation fulfills the following condition: the product of all edge weights around 
a face must equal —1 if the number of edges around the face is mod 4. If the number 
of edges equals 2 mod 4, the product must equal 1 |6] . One can choose an orientation by 
consistently assigning arrows to the edges of the graph, as originally suggested by Kasteleyn 
||TJ|2). The above treatment can be straightforwardly generalized to any genus g Riemann 
surface. 

In the following, we will restrict ourselves to regular 2 x n graphs 02,n embedded on a 
torus, to which we will refer in the following as a strip. On the torus, there are two non- 
trivial cycles, which we will denote by z and w. 

In the case of the plane graph, the edge weights originated solely from the Kasteleyn orien- 
tation. We choose a positive direction on the dimers, say • — > o. Now we assign the weight 
z (w) to each edge which crosses the cycle z (w) in positive direction and the weight \/z 
(1 / w) to each edge which crosses it in negative direction. While the Pfaffian of the Kasteleyn 
matrix yielded a number in the case of the plane graph, it becomes a polynomial in z and w 
on the torus, the so-called characteristic polynomial or Newton polynomial of the graph. The 
coefficient of each monomial z^w^ gives the number of matchings with weight (z, w) = (p, q). 
These are matchings with the number of dimers crossing z in positive direction minus the 
number of dimers crossing z in negative direction equal to p (analogous for q). In the litera- 
ture, what we call the weight is usually referred to as the slope of a height function defined 
on the composition of two matchings^] The matching shown in Figure [l]has weight (1,0), 
where 1 = 2 — 1. 

The partition function, or Newton polynomial, of the dimer model on the torus takes the 
form 

V m , n (z,w) = VfetK = £ N nz ,« a (-iy^+ n ^ Z ^w n \ (1) 

n 2 ,n w 

where the N„ Z/ „ W count the number of matchings of weight (height change) (n z ,n w ). Further- 
more, the total number of matchings for the square graph on the torus is given by 

Z m ,n = \ (-P^(U) + V%k(l, -1) + V$n(-l, 1) + V%k(-1, -1)) , (2) 

1 The height function is defined as follows. Choose a reference matching PMq. To find the slope of a matching 
PM, compose it with the reference matching, PM — PMq, where the minus serves to change the orientation of 
PMq to o — > •. This results in closed loops (composition cycles) and double line dimers. The rule is that when an 
edge in PM belonging to a closed loop is crossed such that the black node is to its left (right), the height changes 
by +1 (—1). If an edge belonging to PMq is crossed, the signs are reversed. This height function is defined up 
to the choice of the reference matching PMq. Crossing the boundary of the fundamental region of the torus, this 
function can jump. If the height function jumps by p units crossing z, it is associated to the power z'' in the 
Newton polynomial of the graph (and equivalently for w). Choosing a different reference matching results in a 
common prefactor of zP°w1° for all monomials. Our method of assigning weights to a matching corresponds to 
choosing a reference matching of weight (0, 0) that does not intersect the z or w cycle. 
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Figure 1: Example of a square graph on the torus 
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n/2 + 1 n/2+x-l n/2+x 
Figure 2: Long strip of squares on the torus 



where the first term is always zero. 



2 Combinatorics for the strip of squares 

We consider a long strip of squares ^2,n on the torus, see Figure |2j focusing on the z-weights. 
A strip containing n black nodes can accommodate matchings with weights n/2, n/2 — 
1, . . . , — (n/2 — 1), —n/2, i.e. there are n + 1 subsets. We would like to find a direct way 
of obtaining the multiplicities of the matchings of a given weight, i.e. the information con- 
tained in the characteristic polynomial, with just the number n as input data. Using the 
Kasteleyn construction, we find the following multiplicities for the first five cases: 
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(3) 



The above sequences do not have an obvious structure. We will solve the problem using an 
operator perspective. 
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(a)d fl (x) (b) di(x,»/2 + *-l) (c) d b (x,n/2 + x) 

Figure 3: The three basic operations that can be performed on a spin pointing up 



2.1 Potts picture derivation 

In order to find a convenient formalism, we use the following picture: we attach to each 
black node a Z m spin, which can point along all the directions in which the black node is 
joined to a white node by an edge. This results in a description reminiscent of the m-state 
Potts model, where m here is the valency of the nodes. A hexagon graph results in a 3-state 
model, while the square graph gives a 4-state model, in which the spin can point up, down, 
left, or right: ||), | j), |— >), |<— ). To describe a dimer configuration, we take the spins to point 
to those white nodes which are joined to the black nodes by a dimer. Since we are interested 
in perfect matchings, we must restrict the possible configurations of the Potts model to those, 
in which each node is only touched by one dimer. 

We label the black nodes by numbers as shown in Figure [Tj The perfect matching shown 
in Figure [Tjcan be written as the following spin state: 



-T-UT) • 



(4) 



As we will show in the following, it is possible to define basic operations on the spins 
allowing us to reach all states starting from the highest weight state. 

The highest weight state is unique and in the spin picture is the one with all spins up: 



ITTT ... TT) =: I"). 



(5) 



Also the lowest weight state is unique and is the one with all spins down. We are now 
looking for basic operations on the spins, which take us from the highest weight state to 
another state. A spin pointing up can be flipped to three possibilities: down, left and right. 
Given the structure of the highest weight matching, flipping a spin down gives another 
perfect matching, while flipping the spin on node X left requires us to flip the spin on node 
n/2 + x — 1 right to arrive at a perfect matching. Flipping the spin on node X right requires 
us to flip the spin at n/2 + x to the left. See Figure |3]for these three basic operations. Each of 
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Figure 4: Long strip of hexagons on the torus 



these operations lowers the weight of the perfect matching by one. We shall denote the spin 
flip from up to down at node x by d a (x). The spin flips by n/2 on the nodes x and y will 
be denoted by d{,(x,y). For x = 2,.. . ,n/2, y = n/2 + x — 1 or y = n/2 + x and for x = 1, 
y = n/2 + x or y = n. 

To familiarize ourselves with this operator picture, see the example of n = 4 black nodes 
in Appendix [A] 

Before we treat the square graph, we first solve the case of the strip of hexagons, which 
can be obtained from the square by deleting half of the vertical links, see Figure |4] 

2.1.1 Combinatorics for the strip of hexagons 

Since the strip of hexagons is a 3-state model and does not accommodate downwards (or 
upwards, depending on which half of the vertical links we choose to delete) pointing spins, 
the weights in z only run from to n /2 and only the spin flips by n/ 2, dj, (x, y), exist. 
The first five cases have the following multiplicities: 

1 2 3 4 5 

n z u z z z° z z 



2 2 1 

4 2 4 1 (6) 
6 2 9 6 1 

8 2 16 20 8 1 

10 2 25 50 35 10 1 

The combinatorics therefore stems exclusively from repeated operations of dj, on the state 
\n). On the hexagon graph, the example of n = 4 black nodes shown in the Appendix 
reduces to what is represented in Figure [5] 

Theorem 2.1. The general formula for the number of perfect matchings with z-weight n/2 — p on 
the strip of hexagons on the torus is 



a 



n,p — . 



F- ,=i ' n-p\ p 
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*IT-T<-> 

Figure 5: Strip of hexagons with n = 4 black nodes 



Proof. The operator description explained above is convenient because it allows to map the 
two-dimensional dimer problem on the hexagonal strip to a one-dimensional monomer- 
dimer problem. In fact we can simply concentrate on the upper line of the strip and consider 
two nodes to be occupied by a dimer if they are occupied by a horizontal line in the strip 
resulting from a b-move, and to be occupied by a monomer if the dimer is vertical in the 
strip. The occupation of the lower line is uniquely determined by consistency. Let us now 
write down the partition function for a monomer-dimer system living on a ID lattice of 
n nodes with periodic boundary conditions (note that we do not require the lattice to be 
bipartite anymore, so n can be odd). The partition function reads: 

L»/2J 

Qn{q) = Yj a n,V°l V > (8) 

where a„ iP is the number of configurations with p dimers. Let us first start with a slightly 
simpler system of n nodes on a line with free boundary conditions, and let P„ (q) be the corre- 
sponding monomer-dimer partition function. This function satisfies the following recursion 
relation: 

Pn+i{q) = Pn{q) + qP„-i{q). (9) 

This can be understood as follows. When adding an extra point after the n-th, one can either 
add a monomer which leaves P n+1 (q) = P n (q) (since P n (q) counts the dimers), or add a 
dimer (and multiply by q) if the last point was previously occupied by a monomer. The 
configurations in P n (q) where the n-th point is a monomer are precisely counted by P n -\ (</)■ 
Note that this is the ^-analogue of the Fibonacci sequence and in fact the actual Fibonacci 
sequence can be recovered by { P„(l) }, where the P n (q) are obtained by using the initial 
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conditions P\(q) = 1, Pi{q) = 1 + q. 

In a similar way, we can understand the relation between P n (q) and Q n (q)- In fact, mak- 
ing the plane strip periodic by adding a line between the first and the n-th node can either 
leave the partition function invariant or add an extra dimer for each configuration where the 
first and last nodes are occupied by a monomer. Therefore, 

Qn(q) = Pn(q) + q Pn-2(q) ■ (10) 

It is easy to see that Q n {q) satisfies the same recursion relation as P n {q)- On the other hand, 
we are interested in the n even case, so it is better to recast it in the form 

Qn+l{q) = (1 + 2q) Q n (q) - q 2 Q„- 2 (q) , (11) 

that can be easily solved with the initial conditions 

Qo(?) = l, Q 2 (q) = l + 2q, (12) 

the solution being 

Qn{q) = 4 K 1 - v^+^)" + f 1 + V^+^Y] ■ (13) 



2" 

To expand p3| in powers of q we can use the identity 

[m/2\ 

E 

k=0 

and find 



[m/2\ / \ 

(a + b) m + {a- b) m = 2 £ J « rf , 



[n/2\ [n/2\ 

E (^) p E 

p=0 A:=p 

Using the identity (see 11201 ) 



(14) 



Q.W - 1 E (;,) i^f = ^ E ( " ) E Q) <««)' • as) 

Inverting the sums, we find 

o.W-^E(«»)'E( 2 " t )C)- 
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we obtain 

Qn(q) = l L^( n Z P )f- (18) 



The explicit expression for Q„ (q) allows to compactly summarize the result for any strip 
by defining the generating function 

00 q — ? 

?^s, q )=Y:s"Q n ( q )= (19) 
Corollary 2.2. The partition function for the strip of hexagons on the torus, considering only z is: 

p hex (z) =z „/2 £(_l)P z -P flB , p = Z " /2 Q„(-1/Z) = 
p=0 



z n/2 
2« 



(20) 



Since there is only one matching with weight it; and one with weight 1/w, the full New- 
ton polynomial on the torus is 

Vg?(z,w)=P* e *(z)-w-±. (21) 

Again we can summarize the result into a generating function in z for the dimers on a 
hexagonal strip as follows: 

2.1.2 Generalization to the square strip 

After having solved the problem for the strip of hexagons, we return to the square strip. We 
denote a state which only contains up and down spins by \a), while a state which contains 
at least one left/right pair by \b). Under the action of our two operators, we have 

d a \a) = \a), d a \b) = \b), d h \a) = \b), d h \b) = \b). (23) 

We denote by p the number of times we act on the highest weight state. The resulting states 
have weight n/2 — p. The order in which the state was acted on by the d a and dy does not 
matter, so at level p, there are p + 1 combinations of d a and dy. On the highest weight state 
\n), we can act in n ways with d a and in n ways with d^. On d a \n), we can then act in n — 1 
ways with d a and in n — 2 ways with dj,. On d(, \ n), we can act in n — 2 ways with d a and in 
n — 3 ways by df,, etc. It is in general easier to compute the number of possibilities of acting 
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with d a on a given state than with dj,, so we choose the most convenient path to obtain the 
full combinatorics, see Figure [6] 
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(n-4)(n-5) \, 
(h-3) X 
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Figure 6: Diagram of d fl/ d;, operations 



We summarize the results for the first three levels: 



V 
1 



d a 

n 
d a d a 

n(n—l) 
2 

da dad a 



+ 



d b 
n 
d a db 



+ / « + 
n (n — 2) 



n(n-l)(w-2) + n(«-2)(n-3) + 



d a d a db 

2! 



dbdb 

n(n-3) 
2 

dbdbda 

n(n— 3)(n— 4) 
2! 



+ 



(24) 



dbdbdb 

n(n— 4)(n— 5) 
3! 



Note the denominators, which avoid an overcounting of states obtained by repeatedly ap- 
plying the same operators. The result for the p-th level is a sum of p + 1 terms. The denom- 
inator of the g-th term on level p is (/! (p — q)\, q = 0, .., p. 
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Theorem 2.3. The q-th term on level p equals to 

1 (n — q — 1)1 n ( p\ ( n — q\ , , 

M IKP-IV- (n-q-p)\ n-q\qj\ p J 

The total number of perfect matchings with weight n/2 — pis therefore 

P t n{n — q — 1)' 
K ' V = k q l(p-q)!(n-q-p)l ' (26) 

Note that & n ,n-p = b n „, which mirrors the symmetry of the sequence. 



Proof We will prove (25 1 by induction. The first three terms are shown explicitly in ( |24| . 
Note that (p — q) counts the number of times d„ has been applied and q counts the number 
of times d;, has been applied. One can act with d a on a given configuration at level {p,q) in 
n — q — p ways, since of the n upwards pointing arrows in the highest weight state, (p — q) 
were turned down by acting with d fl , whereas 2 q were turned horizontally by acting q times 
with d(,. This number has to be normalized by the degeneracy factor q\{p — q)\. Indeed, 

q\(p + 1 - q)\b n/V+ i /q = (n-p-q)q\{p- q)\ b n , Pl q . (27) 

Since every possible state can be arrived at by acting with d a on a state created by acting with 
only db on | n), it is enough to now prove ( |25) for the case p = q, which corresponds to the 
case of the strip of hexagons. In fact b n ,p,p = a, lrP which was calculated in the last section. □ 

Corollary 2.4. The Newton polynomial in zfor the strip with n black nodes (n even) is 

V *\z) = z-"/ 2 t(-l)PzPb n , p = z-" /2 E(-l)PzP E ,, n{n -yV [ (28) 

This constitutes a compact sum formula for the Newton polynomial for the strip of 
squares of arbitrary length. Since there is only one matching with weight w and one with 
weight 1 / w, the full Newton polynomial on the torus is 

<pSfa /W )= V ^ z )- w -L. (29) 

2.1.3 Recursion relation 

Even though we already have the result for the partition function in the form of a sum, we 
will now derive a recursion relation which expresses Vn^(z) through V^_ 1 [z), i.e. we take a 
strip </2,n-2 and add in another square, resulting in Qi iU . Like in the case of the hexagon strip, 
this will provide us with a closed form, and like for the hexagon, the derivation invokes the 
map to a monomer-dimer system. 
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Theorem 2.5. The recursion relation for the partition function of the strip of squares on the torus is 

V:\z) = V 2 s \z) V s n \{z) - V^z) (30) 



with the initial conditions 



p 2 sq (z) =4-z 



(31) 



The solution to ( 30 > is given by 
1 



(-4z) 



n/1 



Z- 1 



l+z(z-6) + z-l + Wl + z(z-6) 



(32) 



This result is equivalent to (28 * and provides a closed form 



Proof. The dimer model on the square strip can be mapped to a one-dimensional monomer- 
dimer in the following way (see Fig. [7|). A dimer living on a z link corresponds to a u- 
monomer living on the • node, a dimer living on a 1 /z link corresponds to a D-monomer 
on the o node, a dimer living on an internal vertical link corresponds to an empty point, 
and a horizontal dimer is again a dimer. In this way, the square strip problem becomes 
a monomer 3 -dimer on a one-dimensional n-node lattice. As for the case of the hexagon 
graph, let us start with free boundary conditions described by the partition function 



P n (u,v,t) = J2 c i,j,ku'v ] t k , 

i,j,k 



(33) 



where c$ is the number of configurations with i u-monomers, j ^-monomers, ft dimers and 
(n — i—j — 2ft) free nodes. 

We can always suppose without loss of generality that the first node is a •, such that 
Pi (u, v, t) = 1 + u. Adding a • node to a string of 2 ft: points results in 



Pik+i(u,v,t) = (l + u)P 2k (u,v,t) +tP 2k -i(u,v,t) , 



(34) 



which can be read as stating that the new point is either free, occupied by a u-monomer or, 
if the 2 ft-th node was free, by a new dimer. Similarly, adding a o node to a string of 2 ft: + 1 
gives 

Pik+i{u,v,t) = (l + v)P 2k+ i{u,v,t) +tP 2k {u,v,t) . (35) 

Adding periodic boundary conditions means adding a new line between nodes 1 and n. One 
can see that the corresponding partition function is given by 



Q n {ll,V,t) = P n (ll,V,t) +tP n - 2 (u,V,t), 



(36) 
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Figure 7: Mapping of a square strip configuration to a monomer-dimer configuration. This 
particular configuration contributes as uv 2 t to the partition function. 



which can be understood as saying that we get either the same configurations or a new dimer 
if the extremal nodes were both empty. Since the relation between Q and P is linear, they 
satisfy the same recurrence equation. On the other hand, Q n is well defined only for n even, 
so the equation is better cast into the form 



Q„ + 2(u,v,t) = [(1 + u) (l + v) + 2t] Q n (u,v,t) - t 2 Q n _ 2 {u,v,t), 



(37) 



and solved with the initial conditions 



Q (u,v,t) = 1, 



Qi(u,v,t) = 1 + u + v + uv + t. 



(38) 



The partition function for the dimer on the square strip is then obtained using the weights 
as they were defined on the initial bipartite graph and reads 



Vn\z) = Q n (-2,-1/2,1). 



(39) 



Hence it satisfies 



KU Z )= (4-^-^n q (z)-^ Sq 2 (2), 



with the initial conditions (pit. Solved explicitly, the closed form (32 1 is obtained. 



(40) 



□ 



The expression in ( [32] ) can be summarized by introducing a generating function in 2 as 
follows. 

4zs (2-1)^2 + 82 



(41) 
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2.2 Comparison with Kasteleyn's product formula 



For the total number of matchings on an m x n square lattice on the torus, a product formula 
exists [2j: 



m/2 n 



£ ^y+b+ab Yl H2\ x 2 sm 2 



a,b=0 



k=l 1=1 



(2k -a) n 



m 



+ x' 2 sin 2 



2l-b)n 



(42) 



Here, x, x' are the edge weights which in our case are equal to one. When specializing 
formula |42]| to m = 2, this gives 



1 " 



Z=l 



2/rZ 



sm 



1 " 



z=i 



(21-1)71 

sm — 



+ 



1 * / . 2/r/ 1 " / . 2 (2/-l)/r 

+^ n 2 v 1 + sm ir + 2 n 2 v 1 + sm 



Z=l 



(43) 



Note, that the first term, A\, is identical to zero (already for general m). Furthermore, Ai = 2, 
and 

A 2 -A t = 2, A 4 -A 3 = 2 Vn. (44) 



From the product formula ( [43| |, this is little obvious, but clear from the point of view of 
formula ^ and the Newton polynomial V^- Since the monomials in w are — w and —1/w, 
^2k ( ' ' 1 ) an< ^ ^zn ( ' ' — ) mus t differ by four. Therefore, we arrive at the identifications 



A 1 = \Vg{\,\), A 2 =\V^ n {l-l), A 3 = \V^ n {-l,l), A 4 = i^(-1 / -1). (45) 



1 T> S 1 



1 ID^l 



1 T)Sq, 



Substituting (45 1 into ( |28| , this means at the same time that 

n(n — q — 1)\ 



2. 



(46) 



p=0 



q=0 



2.3 General case 



We now consider the general case of a square lattice on the torus with n x m nodes (n/2 x 
m/2 black nodes). In general, the states with a given weight have here a much bigger degen- 
eracy, since the interior of the graph allows for many different configurations which do not 
affect the boundaries. Also here, the highest weight state is unique and in the spin picture 
is the one with all spins up. For large examples, it might seem at first surprising that the 
boundary conditions are strong enough to completely fix the configuration in the interior, 
but it is easy to convince oneself by inspecting a small example that the interlacing structure 
of the spins does not allow for any other configuration in the interior. This changes once 
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one of the spins on the boundary points down. For the next simplest case of m = 4 for ex- 
ample, each boundary state with weight n max — 1 has a degeneracy of 3 in the interior. It is 
obvious that this degeneracy grows quickly with growing m and lower weight. Contrary to 
the case before, one must always change more than one spin to arrive at a new matching. 
On the whole we see that the spin picture is not very well adapted to the general case and 
the complexity of the derivation used for the strip rises to an unreasonable level for p > 1. 
Also ansatze for recursion relations have turned out to be of little use. To solve the general 
problem, other methods might be more appropriate. 

3 Translation to the minimal feedback arc set problem 

The problem of finding a minimal set of arcs in a directed graph upon the deletion of which 
the graph becomes acyclic is well studied in mathematics and computer science [l2"T| . Such 
a set is called a minimal feedback arc set (FAS). There exists a precise relation between the 
problem of finding all minimal FAS of a digraph and the problem of identifying all the perfect 
matchings in its dual graph. 

Consider a bipartite plane graph 'S with N nodes, and its graph dual < S' (where nodes 
become faces, faces become nodes and edges remain edges) . The dual graph <$' becomes a 
digraph if the edges around a face corresponding to a black node in are oriented clockwise, 
while the edges circling a face corresponding to a white node are oriented counterclockwise. 
The one-cycles in e S' are generated by the plaquettes { pj which correspond to the ver- 
tices of . Removing the edge e,y shared by the cycles p,- and pj breaks both cycles. A minimal 
FAS is obtained taking the collection of N/2 edges { e^ k },? shared by disjoint pairs of pla- 
quettes p{. and pi k . In the dual graph <S , this corresponds to selecting a set of edges joining 
all the nodes and touching all of them only once. In other words, a minimal FAS in is a 
perfect matching in <$ . 

The situation is different when # is embedded on a Riemann surface of genus g > 0, 
because in addition to the one-cycles generated by the plaquettes, there are 2g equivalence 
classes of cycles of non-trivial holonomy. In this case, the winding cycles in ( S' are gener- 
ated by the zig-zag paths^] It follows that being a perfect matching in <& is only a necessary 
condition for a set of edges to be a FAS in ( S' . Let us again restrict ourselves to the case of 
g = 1. A useful way to represent the partition function in Eq. ([T]l consists in drawing a point 
in the (z, w) plane at coordinates (n z , n w ) for each monomial of the form z" z w nw , correspond- 
ing to the matchings with weight (n z , n w ). In this way, one obtains the region of a Z 2 lattice 
delimited by a convex polygon, the Newton polygon. We can hence distinguish between in- 
ternal and boundary points (or matchings). It was shown in |Tl8l that a perfect matching 
corresponding to an internal point of this Newton polygon is always a FAS. Removing a 

2 A zig-zag path is a path which turns alternatingly maximally left and maximally right at the vertices. 
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Figure 8: Feedback arc sets on the dual graph to the square strip. The dual graph is rep- 
resented in light blue. The removal of the boundary matchings (dashed lines) preserves a 
zig-zag path (dark blue). 



boundary matching, on the other hand, always preserves at least one zig-zag path. 

In the case at hand, i.e. for the square strip, the polygon is a rhombus with the four 
corner vertices at the points ( ±n / 2, 0) and (0, ± 1 ) . This means that there are four boundary 
matchings that will not be feedback arc sets. In fact, removing the arrows corresponding to 
these matchings will always preserve exactly one zig-zag path (see Fig.lBl. 



Theorem 3.1. On a 2 x n digraph on the torus (n even), there are 



1 + V2) +I-1-V2 



minimal feedback arc sets. This number is generated by the function 

2 ( 1 + y/2] + ( \fl - l) s 

F^( S ) = £n„ fa v = 

n=0 



(47) 



1) (-1+ (-I + V2 



r^s,-i)- Y 



(48) 



Proof. Follows from combining Eq. (|2j and Eq. ( p2| and subtracting the four boundary match- 
ings. □ 
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4 Discussion and further directions 



In this note, we have derived an explicit formula for the partition function of the dimer 
model on a 2 x n strip of squares or hexagons on the torus. This result was derived via an 
operator picture, and by mapping the problem to a one-dimensional monomer-dimer sys- 
tem. This result was furthermore translated to the question of finding the minimal feedback 
arc sets on the dual graph. 

An obvious continuation of this work would be the generalization to an n x m graph. 
This problem turns out to be much more complicated than the one addressed here and the 
methods employed here have proven not to be properly adapted to the more general ques- 
tion. 
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A Example: n = 4 square strip 



To familiarize ourselves with this operator picture, we consider the example of n = 4 black 
nodes, see Figure [9] We start with the highest weight state, which has has weight 2. The 
solid arrows in the diagram denote the action by d fl , the dotted arrows the action by dj,. We 
find a symmetric sequence structure. 



mm | 




'%i if 



Figure 9: Strip with n = 4 black nodes 
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